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High  concentration  photovoltaic  technology  promises  the  large-scale  generation  of  clean-renewable 
energy  with  competitive  costs.  Like  any  other  systems  for  electricity  generation,  it  is  important  to  know 
the  electrical  characteristics  of  the  system.  However,  while  there  is  a  wide  experience  in  modeling  the 
behavior  of  traditional  photovoltaic  systems,  not  every  model  for  flat-plate  solar  cells  or  modules  is 
directly  applicable  to  high  concentration  photovoltaic  cells  or  modules  because  of  the  special  features  of 
these  devices  (use  of  multijunction  cells,  use  of  optics  for  high  concentration,  etc.).  So,  in  recent  years,  the 
scientific  community  has  devoted  considerable  efforts  in  developing  models  that  reproduce  the  electrical 
behavior  of  high  concentration  cells  and  modules.  These  models  allow  calculating  the  main  electrical 
parameters  of  the  device  from  its  operating  conditions  (irradiance,  cell  temperature,  spectral  distribution 
of  the  radiation,  etc.).  In  this  paper,  a  comprehensive  review  of  existing  models  for  the  electrical 
characterization  of  high  concentration  photovoltaic  cells  and  modules  is  presented  with  the  aim  of 
helping  the  photovoltaic  professionals  and  researchers  in  the  design,  monitoring  and  energy 
prediction  tasks. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

High  Concentration  Photovoltaic  (HCPV)  cells  and  modules  oper¬ 
ate  under  concentration  ratios  between  300  and  2000  suns.  The 
concentration  mechanism  is  done  by  lenses  or  mirrors  that  either 
reflect  or  refract  the  sunlight  on  the  solar  cell  surface.  Replacing  the 
expensive  solar  cell  materials  by  cheaper  optics  offers  the  opportu¬ 
nity  of  reducing  the  costs  of  the  modules.  In  addition,  the  use  of 
multijunction  solar  cells  implies  an  increase  of  the  system  efficiency, 
what  also  can  reduce  the  associated  costs  [1], 

Three  main  elements  are  required  in  order  to  generate  elec¬ 
tricity:  the  multijunction  solar  cell,  the  optics  and  the  solar  tracker. 
These  elements  can  be  combined  in  a  lot  of  different  configura¬ 
tions  [2,3],  A  typical  commercial  HCPV  system  is  composed  of 
concentration  modules  mounted  on  two-axis  solar  trackers,  where 
each  module  contains  several  multijunction  solar  cells  and  each 
cell  has  its  own  optics,  typically  a  point-focus  Fresnel  lens.  Triple¬ 
junction  GalnP/GalnAs/Ge  cells  are  usually  utilized  [4],  while  other 
material  combinations  are  now  under  research  [5-8].  Other 
auxiliary  elements  are  used,  such  as  heat  sinks,  Secondary  Optical 
Elements  (SOEs),  bypass  diodes  and  protection  elements. 
Nowadays,  these  configurations  dominate  the  HCPV  market  and 
are  very  promising  solutions  for  the  future  of  photovoltaics  [9], 

Modeling  the  electrical  behavior  of  HCPV  cells  and  modules  is 
an  important  task  for  the  design,  monitoring  and  energy  predic¬ 
tion  of  these  kinds  of  systems.  In  this  paper,  a  review  of  the  state- 
of-the-art  of  models  that  reproduce  the  operation  of  HCPV  ceils 
and  modules  is  presented.  These  models  calculate  the  l-V  char¬ 
acteristic  curve  of  the  device  from  its  operating  conditions.  Not 
every  model  allows  calculating  the  whole  I-V  curve;  some  of  them 
only  calculate  some  relevant  electrical  parameters,  such  as  the 
maximum  power  (PMpp),  the  short-circuit  current  (/sc),  the  open- 
circuit  voltage  (Voc),  the  fill  factor  (FF),  etc. 

The  presented  models  differ  in  the  parameters  required  to 
characterize  the  device,  in  the  variables  used  as  inputs  and,  finally, 
in  the  equations  and  procedures  used  to  obtain  the  electrical 
characteristics.  The  models  present  different  levels  of  complexity 
and  accuracy  and  require  different  equipments  to  be  applied.  This 
review  remarks  the  characteristics  of  each  one  of  these  models  in 
order  to  help  in  the  choice  of  the  most  suitable  model  for  each 
specific  application. 

The  paper  is  organized  as  follows:  in  the  second  section,  a  brief 
overview  of  the  technologies  under  study  is  presented;  in  the  third 
section,  the  existing  models  for  multijunction  solar  cells  are 
summarized;  the  fourth  section  describes  the  main  contributions 
in  the  modeling  of  HCPV  modules;  finally,  section  five  comments 
the  conclusions  of  the  work. 


2.  Brief  technology  overview 

The  main  elements  that  constitute  a  HCPV  generator  are  the 
solar  cells,  the  receivers,  the  optical  devices  and  the  tracking 
system.  A  collection  of  solar  receivers  with  electrically  intercon¬ 
nected  cells  and  with  integrated  optical  devices  is  known  as  a 
HCPV  module.  Next,  the  main  characteristics  of  these  elements  are 
described. 


2.3.  Solar  cells 

High  concentration  solar  cells  have  specific  characteristics 
different  to  those  of  conventional  flat-plate  solar  cells.  The  HCPV 
cells  are  designed  to  extract  more  energy  under  a  larger  solar 
radiation.  These  cells  have  to  be  exclusively  made  of  high  quality 
semiconductor  materials  because  of  the  significance  of  the  con¬ 
version  efficiency,  a  factor  which  is  directly  related  to  the  quality 
of  the  crystal  used,  in  which  impurities  must  be  avoided. 

An  important  characteristic  of  these  solar  cells  is  their 
enhanced  thermal  behavior.  Under  high  temperatures,  both  Si 
and  GaAs  cells  suffer  a  significant  decrease  of  performance  mainly 
because  of  the  voltage  reduction  produced  by  a  temperature 
increase.  Nonetheless,  this  effect  is  less  accurate  at  high  concen¬ 
tration  levels. 

Despite  the  fact  that  many  different  types  of  cells  are  being 
tested  at  laboratory  level  and  although  the  reported  efficiency  has 
increased  over  40%  for  multijunction  cells  [10],  many  of  the 
concentration  systems  installed  nowadays  use  silicon  solar  cells, 
providing  efficiencies  which  hardly  exceed  25%  under  concentra¬ 
tion.  However,  most  of  the  current  developments  are  based  on 
triple-junction  cells.  These  cells  can  be  described  as  a  stack  of  cells 
of  different  compositions  with  a  progressively  decreasing  band- 
gap,  so  that  the  cell  on  the  top  of  the  stack  absorbs  higher  energy 
photons  than  the  cell  in  the  middle,  and  this  absorbs  higher 
energy  photons  than  the  cell  in  the  bottom. 

2.2.  Receivers 

Solar  cells  have  to  be  mounted  on  a  receiver  to  extract  the 
generated  current  and  to  dissipate  and  remove  the  heat  resulting 
from  sun  power.  Current  extraction  is  made  by  thick  wires  or 
ribbons  that  are  soldered  to  the  cells,  whereas  heat  removal  is 
carried  out  by  a  heat  sink  added  to  the  receiver.  Materials  like 
copper  or  aluminum  are  currently  used  for  heat  evacuation.  Also, 
the  solar  receivers  often  include  a  bypass  diode  connected  in 
parallel  with  the  cell  in  order  to  avoid  overheating  of  the  cell  in 
non-usual  operating  conditions  such  as  shading. 

2.3.  Optical  devices 

The  aim  of  the  optical  devices  used  in  this  technology  is  to 
concentrate  radiation  and  to  increase  the  luminous  flux  on  the 
cells.  Optical  devices  may  be  simple  or  consisting  of  primary  and 
secondary  optical  element  (SOE).  Primary  optical  devices  collect 
direct  sun  rays.  Secondary  optical  devices  receive  the  light  from 
the  primary  ones  and  then  change  its  properties  by  means  of 
spectral  filtering,  light  homogenization  or  light  direction  changing. 
In  this  respect,  different  optical  devices  are  currently  used:  circular 
parabolic  dishes,  parabolic  dishes  with  secondary,  square  flat 
Fresnel  lenses,  square  flat  Fresnel  lenses  with  secondary,  linear 
flat  lenses,  linear  arched  lenses  or  linear  parabolic  reflectors. 

Most  of  the  systems  currently  designed  are  based  on  Fresnel 
lenses.  This  kind  of  lens  reduces  the  amount  of  material  required  to 
concentrate  the  light  by  splitting  the  lens  into  a  set  of  concentric 
annular  sections  known  as  Fresnel  zones.  The  use  of  these  zones 
allows  keeping  the  required  curvature  without  increasing  the 
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thickness  by  means  of  adding  discontinuities  between  them.  Fig.  1 
shows  an  example  of  a  solar  receiver  integrated  with  point-focus 
Fresnel  lens. 

2.4.  HCPV  module 

A  HCPV  module  is  the  smallest,  complete,  environmentally 
protected  assembly  of  receivers  and  optical  devices  that  is  able  to 
transform  an  input  of  unconcentrated  solar  radiation.  Consequently, 
a  typical  HCPV  module  is  made  up  of  a  group  of  cells,  primary  optics, 
secondary  optics  (optional)  and  housing  components  such  as 


interconnection  and  mounting.  Fig.  2  shows  an  example  of  a  typical 
HCPV  module  and  its  components. 

2.5.  Tracking  system 

Point-focus  based  HCPV  modules  must  be  always  mounted  on 
two-axis  solar  trackers,  i.e.  the  modules  must  be  always  pointing 
the  solar  rays  in  order  for  the  lenses  to  be  able  to  focus  the 
radiation  on  the  small  solar  cell  area.  This  kind  of  tracking  system 
is  more  complex  than  linear  tracking  or  static  systems  from  a 
mechanical  point  of  view  but  it  presents  the  advantage  of  max¬ 
imizing  the  solar  radiation  capture  along  the  daily  trajectory. 


Primary 

concentrator 


Solar 

receiver 


3.  HCPV  cells  models 

In  this  section,  five  existing  models  that  reproduce  the  behavior 
of  a  multijunction  solar  cell  operating  under  high  concentration 
are  summarized.  The  presented  models  assume  certain  simplifica¬ 
tions  in  order  to  be  applicable  in  engineering  studies.  Other  more 
sophisticated  models  that  divide  the  cell  into  hundreds  or  thou¬ 
sands  of  elementary  cells  are  not  included  in  this  review  [11,12]. 
These  distributed  models  are  difficult  to  implement  and  require 
high  computational  effort,  so  that  they  are  not  commonly  used  at 
the  engineering  level. 


Fig.  1.  Schematic  of  a  HCPV  solar  receiver  integrated  with  point-focus  Fresnel  lens.  3.1.  Models  based  on  the  equivalent  circuit  of  the  multijunction  cell 


Fig.  2.  Schematic  of  a  HCPV  module  based  on  point-focus  Fresnel  lenses. 
Source:  IEC  62108  [2], 


A  triple  junction  solar  cell  can  be  considered  as  composed  of  three 
series-connected  p-n  junctions.  Each  junction  can  be  represented  by 
an  equivalent  circuit  model,  being  the  most  used  models  the  Single 
Exponential  Model  (SEM)  and  the  Double  Exponential  Model  (DEM). 
Both  models  differ  in  the  number  of  diodes  that  characterize  the 
saturation  current.  By  connecting  the  three  subcells  circuits  in  series, 
the  equivalent  circuit  for  the  triple  junction  cell  is  obtained,  Fig.  3. 

Taking  into  account  the  DEM  model,  the  cell  I-V  curve  can  be 
mathematically  expressed  by  the  next  set  of  equations 

v=  i  vt  (i) 

i=  1 

1  =  fpH,i-foi,i{exp[(Vj  +  IRs,i)/mi,iVT]-l} 

-/02,i{exp[(Vi  +  lRs,i)/m2iVT] 

-l}-(Vi  +  IRSii)/Rp,i  ,  f=l,  2,  3  (2) 

where  the  subscript  i  indicates  the  subcell  (l=top,  2= 
middle,  3  =  bottom).  Each  subcell  is  characterized  by  a  set  of 
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Fig.  3.  Equivalent  circuit  models  for  triple  junction  solar  cells  by  using  the  Single  Exponential  Model  (left)  and  the  Double  Exponential  Model  (right). 
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7  parameters:  /pHl!  (the  photocurrent),  /01i,  and  I02,i  (the  saturation 
currents),  irq ,  and  m2,;  (the  diode  ideality  factors),  Rsi  (the  series 
resistance)  and  RPi  (the  shunt  resistance).  VT  is  the  semiconductor 
thermal  voltage,  Vr=kT/q,  being  k  the  Boltzmann  constant,  T  the 
absolute  temperature  of  the  junction  and  q  the  electron  charge. 
The  equations  for  the  SEM  scheme  are  as  follows: 

V=  |  (3) 

i  =  1 


1  =  fpH,i-fo,i{exp[(V,-  +  VT]-1 } 

-(Vi  +  /Rs,i)/%  ,  i=l,  2,  3  (4) 

where  each  subcell  is  now  characterized  by  5  parameters:  IPHi  (the 
photocurrent),  I0j  (the  saturation  current),  m,-  (the  diode  ideality 
factor),  Rsi  (the  series  resistance)  and  RPi  (the  shunt  resistance). 

In  general,  the  model  parameters  depend  on  the  operating 
conditions  of  the  multijunction  cell,  such  as  concentration  and  cell 
temperature.  However,  as  multijunction  cells  are  often  monolithic 
devices  which  incorporate  only  two  electrical  terminals,  the  sepa¬ 
rate  measurement  of  the  individual  subcells  is  not  possible.  There¬ 
fore,  the  determination  of  the  subcells  parameters  is  very  difficult. 

Some  authors  have  proposed  methods  for  extracting  the  sub¬ 
cells  parameters,  or  at  least  some  of  them,  by  applying  different 
assumptions  and  simplifications.  We  can  find  several  extraction 
methods  based  on  both  the  SEM  scheme  [13-18]  and  the  DEM 
scheme  [17,19-23],  However,  it  is  remarkable  that  these  studies 
often  estimate  parameters  only  at  particular  operating  conditions 
and  not  for  the  whole  range  of  possible  conditions,  so  that  the 
application  of  these  models  is  sometimes  limited. 

3.2.  Model  of  GS.  Kinsey  et  al.  (2008) 

The  model  proposed  by  G.S.  Kinsey  et  al.  [24,25]  allows 
calculating  the  maximum  power  (PMpp)  of  a  triple  junction  solar 
cell  from  its  operating  conditions:  concentration  ratio,  X,  cell 
temperature,  Tc,  and  incident  spectrum,  £'(/).  First,  three  electrical 
parameters  of  the  solar  cell  as  a  whole  are  calculated:  the  open- 
circuit  voltage,  Voc,  the  fill  factor,  FF,  and  the  short-circuit  current 
density,  Jsc.  Once  these  parameters  have  been  determined,  the 
PMpp  is  obtained  by: 

Pmpp  =  FFJsc.4celVoc  (5) 

ACel  being  the  solar  cell  active  area. 

The  Voc  and  FF  are  assumed  to  be  independent  of  the  incident 
spectrum.  While  this  assumption  is  adequate  for  the  Voc  [26,27], 
the  FF  variations  with  spectrum  can  have  a  significant  impact  on 
performance  [28,29].  However,  in  this  model,  these  variations  are 
neglected.  The  V0c  is  obtained  from  the  reference  Voc  at  one  sun 
and  25  °C  of  cell  temperature  by  applying  a  logarithmic  correction 
with  concentration  and  a  linear  correction  with  temperature. 
On  the  other  hand,  the  FF  is  obtained  from  the  reference  FF  by 
only  applying  a  linear  correction  with  temperature. 

The  model  requires  knowing  the  External  Spectral  Response  of 
the  component  subcells  (ESR,)  as  a  function  of  the  wavelength,  X, 
at  different  cell  temperatures.  These  functions  must  be  measured 
with  the  help  of  a  solar  simulator.  A  solar  simulator  usually 
consists  of  three  major  components:  (1)  light  source(s)  and 
associated  power  supply;  (2)  any  optics  and  filters  required  to 
modify  the  output  beam  and  (3)  the  necessary  controls  to  operate 
the  simulator,  adjust  irradiance,  etc.  [30]. 

The  model  equations  for  calculating  the  short-circuit  current 
density  are: 


-/pHmid  —  J  ESRmid (2,  T c  )E(A)dA  (7) 

Jsc  =  mint/pHtop  JpHmidl  (8) 

where  /PHtoP.  ./pHmid  are  the  photocurrents  of  the  top  and  the 
middle  subcells.  Eq.  (8)  considers  that  the  short-circuit  current  of 
the  cell  is  limited  by  the  smaller  photocurrent  of  the  top  and 
middle  subcells,  i.e.  the  bottom  subcell  never  limits  the  current. 

The  model  requires  knowing  several  parameters  that  charac¬ 
terize  the  multijunction  cell,  being  the  most  important  the  ESR  of 
the  top  and  the  middle  subcells.  These  values  are  obtained  with 
the  help  of  the  HIPSS  solar  simulator  from  Spectrolab  [31].  In 
addition,  the  use  of  a  spectroradiometer  is  required  to  measure  the 
incident  spectrum.  Alternatively,  the  spectrum  can  be  estimated 
from  atmospheric  parameters  by  using  simulation  software  such 
as  SMARTS  [32,33], 


3.3.  Model  of  C.  Dominguez  et  al.  (2010) 


The  direct  measurement  of  the  spectrum  is  difficult  and 
requires  the  use  of  a  spectroradiometer,  which  is  not  always 
available.  Another  option  is  to  use  isotype  cells.  These  sensors 
have  the  same  spectral  response  than  the  component  subcells  of 
the  multijunction  device.  This  way,  the  photocurrents  measured  by 
the  isotype  cells  under  the  same  spectrum  than  the  multijunction 
cell  allow  knowing  the  photocurrents  of  the  component  subcells. 

The  model  proposed  by  C.  Dominguez  et  al.  [16]  is  based  on  the 
use  of  isotype  cells.  It  is  a  translation  model,  which  calculates  the 
points  of  the  I-V  curve  at  the  desired  operating  condition  from  the 
points  of  the  I-V  curve  at  the  reference  condition.  Three  transla¬ 
tions  are  applied  to  each  point  of  the  reference  I-V  curve:  (1)  a 
translation  of  the  current  coordinate  to  the  desired  value  of 
concentration:  (2)  a  translation  of  the  voltage  coordinate  to  the 
desired  values  of  concentration  and  incident  spectrum;  (3)  a 
correction  of  the  voltage  coordinate  to  take  into  account  the 
change  in  the  cell  operating  temperature.  This  way,  the  model 
calculates  the  I-V  characteristics  as  a  function  of  irradiance, 
spectrum  and  cell  temperature. 

Translation  equations  are  analytically  obtained  by  considering 
the  SEM  equivalent  circuit  model  for  multijunction  cells  (Fig.  3). 
In  this  scheme,  the  three  shunt  conductances  are  neglected  and 
the  three  series  resistances  are  lumped  into  a  single  effective 
series  resistance,  Rs.  For  instance,  the  translation  of  the  voltage 
coordinate  due  to  the  change  in  concentration  and  spectrum  is 
represented  by  the  following  equation 


V(Tcref)  =  Vref  + 


kjTcref  +  273.15) 


mMjln(X)+ln 


Ftop/PHtop,ref-freAmt0|VFi 
^PHtop,ref— ^ref 


mkMPHmidjef  ^ref 
^PHmid,ref— ^ref 


F  hnt^Pl 


botiPHbot,ref  ^ref 
^PHbot,ref— ^ref 


J-/ref(X-l)Rs 


(9) 


(Vref,  /ref)  being  the  considered  point  of  the  reference  I-V  curve, 
V(TCref)  the  translated  voltage,  /pHi.ref  the  photocurrents  of  the 
component  subcells  at  reference  conditions,  X  the  desired  con¬ 
centration,  m,  the  diode  ideality  factors  of  the  component  subcells 
and  mMj  an  effective  ideality  factor  for  the  multijunction  device. 
The  factors  F,-  are  defined  for  each  subcell  as 

p _ ZpHi //pHi.ref _  (10) 

min{/pHtop,  /pilmid- /pHbot}/ttlin{/pHtop,refi  /pHmid, ref  ■  /pHbot.refl 


./pHtop  —  J  ESRtop(2,  Tc)E(A)dA 


i.e.  they  account  for  the  photocurrents  of  the  subcells,  IPHi  as 
measured  by  the  isotype  cells. 
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The  model  depends  on  some  parameters  that  characterize  the 
multijunction  device.  Some  of  them  can  be  found  in  semiconductor 
textbooks  but  some  other  must  be  fitted  by  using  simulation  tools 
that  minimize  the  error  between  experimental  data  and  modeled 
data.  A  solar  simulator  is  required  to  obtain  the  experimental  data. 


3.4.  Model  of  E.F.  Fernandez  et  al.  (2013) 


The  aim  of  the  model  proposed  by  E.F.  Fernandez  et  al.  [34,35] 
is  to  obtain  a  model  based  on  easily  adjustable  parameters. 

The  measurements  of  several  multijunction  cells  showed  that 
the  Ge  bottom  subcell  never  limits  the  short-circuit  current  of  the 
multijunction  device  and  has  a  small  influence  on  the  rest  of  the 
electrical  parameters  under  normal  operating  conditions.  There¬ 
fore,  the  behavior  of  the  cell  can  be  modeled  by  considering 
only  two  “equivalent  subcells”  that  verify  the  equation  of  the 
SEM  model. 

The  short-circuit  currents  of  the  two  equivalent  subcells  (/Sci. 
hc2)  are  supposed  to  be  equal  than  the  photocurrents  of  the  top 
and  the  middle  real  subcells.  These  photocurrents  must  be 
measured  with  isotype  cells  or  calculated  from  the  incident 
spectrum  as  measured  by  a  spectroradiometer  if  the  ESR  of  the 
subcells  is  known. 

On  the  other  hand,  the  open-circuit  voltages  of  the  equivalent 
cells  (Vocl,  VOC2)  verify  the  condition  that  their  sum  is  equal  to  the 
open-circuit  voltage  of  the  real  cell.  A  logarithmic  correction  on 
concentration  is  used  to  calculate  these  open-circuit  voltages  from 
the  V0c  values  at  one  sun.  The  objective  of  calculating  these  open- 
circuit  voltages  is  to  obtain  the  saturation  currents  of  the  equiva¬ 
lent  cells  (/0i,  I02),  which  are  computed  by  using  well-known 
simplified  expressions  derived  from  the  SEM  solar  cell  model. 

The  equation  of  the  I-V  curve  of  the  multijunction  cell  is 
modeled  by  means  of  a  series  association  of  the  two  equivalent 
subcells,  i.e.  the  voltages  of  these  subcells  are  summed  in  order  to 
obtain: 


V  = 


mk(7c  +  273.15) 
<7 


+  In 


-IRs 


OD 


Rs  and  m  being  effective  values  for  the  series  resistance  and  the 
diode  ideality  factor  of  the  multijunction  device.  Finally,  this  I-V  curve 
is  corrected  taking  into  account  an  effective  shunt  resistance  RP. 

Because  of  the  simplicity  of  the  model  equations,  the  imple¬ 
mentation  of  the  numerical  algorithm  for  obtaining  the  model 
parameters  is  easy.  However,  like  in  the  models  commented  above, 
a  set  of  measurements  carried  out  in  a  solar  simulator  is  required 
for  adjusting  these  parameters. 

3.5.  Model  of  Sandia  National  Laboratories  (2004) 

The  models  commented  above  cannot  always  be  applied  in 
practice,  because  of  the  need  of  a  solar  simulator,  the  need  of 
instruments  such  as  spectroradiometer  or  isotype  cells  or  the 
difficulty  of  implementing  the  algorithms  for  adjusting  the  model 
parameters.  So,  in  some  cases,  the  use  of  a  simpler  model  could  be 
justified  in  spite  of  a  little  loss  in  accuracy.  The  model  developed  at 
Sandia  National  Laboratories  [36]  was  mainly  conceived  for 
characterizing  modules  and  arrays  of  different  technologies,  but 
it  also  can  be  applied  to  high  concentration  solar  cells.  It  does  not 
require  sophisticated  instruments  and  the  model  parameters  can 
be  obtained  from  outdoor  measurements. 

As  it  has  been  shown,  the  main  problem  for  characterizing 
multijunction  cells  is  that  they  are  very  sensitive  to  the  incident 
spectrum.  In  general,  the  sunlight  spectrum  varies  as  a  function  of 
air  mass,  clouds,  turbidity  due  to  aerosol  effects  and  precipitable 
water  vapor  [28].  However,  the  main  contribution  to  the  spectrum 


change  is  air  mass  [37,38].  In  the  model  of  Sandia  National 
Laboratories,  the  measurement  of  the  sunlight  spectrum  is  not 
required  and  the  spectral  influences  are  only  quantified  by  the  air 
mass  parameter,  as  defined  in  [39]. 

The  model  equations  are  formulated  as  a  function  of  the 
“effective  irradiance”  (Bef)-  In  the  case  of  high  concentration  cells, 
this  parameter  is  defined  as: 

Bef  =  [GB/,(AM)]/GS  (12) 

Gb  being  the  incident  direct  irradiance  and  CB*  the  reference 
direct  irradiance.  /i(AM)  is  the  spectral  correction  factor  as  defined 
in  the  standard  ASTM  E  973  [40],  which  is  approximated  by  a 
fourth  order  polynomial  function  of  the  air  mass. 

Five  points  of  the  cell  I-V  curve  are  calculated  by  the  model:  the 
three  classical  points  (short-circuit,  open-circuit  and  maximum 
power  points)  and  two  additional  points:  (l/x.  be).  the  point  at  a 
voltage  Vx=V0c/2,  and  (VxxJxx).  the  point  at  a  voltage  Vxx= 
(V'mpp+V'oc)/2.  Fig.  4  indicates  the  position  of  these  points  on  the 
cell  I-V  curve.  As  an  example  of  the  model  equations,  the 
expressions  that  allow  calculating  the  ISc  and  V0c  parameters  are 
the  following 


he  =  BefUsc  +  «isc(Tc-T*)]  (13) 


Voc  =  VSc  +  <5ln(Bef)  +  A/oc^c-T?)  04) 

where  the  superscript  indicates  parameters  at  reference  con¬ 
ditions.  In  these  equations,  the  influence  of  cell  temperature  is 
modeled  with  the  temperature  coefficients  «|SC  and  /?Voo  where 
Pm oc  is  allowed  to  vary  linearly  with  the  effective  irradiance. 
Coefficient  S  is  the  product  of  the  effective  diode  ideality  factor, 
m,  and  the  thermal  voltage  VT. 

Every  model  parameter  can  be  obtained  by  regression  analysis 
from  outdoor  monitorized  data,  what  is  an  advantage  with  respect 
to  other  models.  The  main  disadvantage  of  this  model  is  that  it 
does  not  take  into  account  ambient  factors  that  affect  the  sunlight 
spectrum,  such  as  clouds,  aerosols  and  precipitable  water  vapor. 
The  model  is  oriented  to  clear  days  with  clean  atmosphere. 
However,  it  is  remarkable  that  these  conditions  are  prevalent  at 
the  locations  where  HCPV  systems  are  installed. 

Table  1  summarizes  the  main  characteristics  of  the  referenced 
models  for  high  concentration  photovoltaic  cells,  including  the 
required  inputs  and  advantages  of  each  one. 


Fig.  4.  Points  of  the  I-V  curve  calculated  by  the  model  of  Sandia  National 
Laboratories. 
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4.  Models  for  estimating  the  maximum  power  of  a  HCPV 
module 

This  section  summarizes  the  state-of-the-art  of  models  for 
HCPV  modules  characterization.  These  models  take  into  account 
different  atmospheric  parameters  that  affect  the  behavior  of  the 
module,  such  as  direct  normal  irradiance  (DNI),  air  mass  (AM),  air 
temperature  (Tair),  wind  speed  (WS),  etc.,  and  some  of  the  models 
use  the  mean  cell  temperature  of  the  module  (Tc)  as  input. 

4.1.  Model  of  the  standard  ASTM  E  2527 

The  standard  ASTM  E  2527  [41  ]  defines  a  very  simple  way  of 
estimating  the  Pmpp  of  a  high  concentration  module.  It  uses  as 
inputs  the  direct  normal  irradiance  (DNI),  the  air  temperature  (Tair) 
and  the  wind  speed  (WS).  It  is  based  on  a  single  equation: 

Pmpp  =  DNI(Ci  +  C2DNI  +  C3Tajr  +  C4WS)  (15) 

where  Cj,  c2,  c3  and  c4  coefficients  are  obtained  by  regression 
analysis  of  outdoor  monitorized  data.  As  can  be  seen,  the  model 
does  not  take  into  account  the  spectral  effects  that  significantly 
affect  the  behavior  of  modules  containing  multijunction  solar  cells. 

4.2.  Model  of  G.  Peharz  et  al.  (2011 ) 

The  model  of  G.  Peharz  et  al.  [42]  does  take  into  account  the 
spectral  influences.  It  is  based  on  the  following  equation: 

Pmpp  =  cDNiDNI  +  cZ2Z2  +  czZ  +  cICTc  +  coffset  (16) 

where  the  Z  parameter  represents  the  mismatch  between  the 
photocurrents  of  the  top  and  middle  subcells  with  respect  to  the 
reference  AM1.5d  spectrum.  The  value  of  Z  is  obtained  from  the 
photocurrents  as  measured  by  isotype  cells  [43].  A  value  of  Z>  0 
corresponds  to  a  blue  spectrum,  where  the  limiting  subcell  is  the 
middle  subcell,  while  a  value  of  Z  <  0  corresponds  to  a  red 
spectrum,  where  the  limiting  subcell  is  the  top  subcell.  In  order 
to  quantify  the  value  of  Z,  some  parameters  of  the  solar  cell  must 
be  measured  in  a  solar  simulator.  cDNi,  Cz2,  cz,  cTc  and  c0frset 
coefficients  are  obtained  by  regression  analysis  from  outdoor 
monitorized  data. 

4.3.  Model  of  M.  Steiner  et  al.  (2012) 

The  model  proposed  by  M.  Steiner  et  al.  [44]  is  a  complete 
model  that  combines  the  spectrum  simulation  software  SMARTS 
[32,33]  with  ray  tracing  software  for  modeling  the  optics  of  the 
Fresnel  lenses  [45]  and  a  SPICE  network  model  for  modeling  the 
l-V  characteristics  [46,47].  The  model  requires  a  complete  set  of 
instruments  in  order  to  be  applied  because  it  takes  into  account 


many  factors  that  affect  the  module  behavior.  Beyond  a  pyrheli- 
ometer  and  weather  station  equipment,  it  requires  isotype  cells  in 
order  to  measure  the  Z  parameter,  multifilter  rotating  shadowband 
radiometer  for  the  determination  of  Aerosol  Optical  Depth  (AOD) 
and  Precipitable  Water  (PW)  and  tracking  accuracy  sensor  for 
measuring  the  alignment  of  the  tracker  to  the  sun.  The  model 
allows  calculating  the  Pmpp  of  high  concentration  modules  with  a 
high  level  of  accuracy,  being  its  main  disadvantages  the  difficulty 
of  implementing  and  the  equipment  requirements. 

4.4.  Model  of  A.J.  Rivera  et  al.  (2013) 

Due  to  the  fact  that  the  relation  between  atmospheric  para¬ 
meters  and  module  output  maximum  power  is  complex  and  non¬ 
linear,  the  use  of  artificial  neural  network-based  models  has  been 
proposed.  This  approach  has  been  followed  by  A.J.  Rivera  et  al.  [48]. 
In  particular,  a  cooperative-competitive  hybrid  algorithm  for  radial 
basis  function  networks  was  implemented. 

The  model  uses  the  Average  Photon  Energy  (APE)  [49,50]  in 
order  to  quantify  the  spectral  influences  on  HCPV  module  max¬ 
imum  power.  This  is  a  single  value  in  eV  that  characterizes  the 
shape  of  the  spectrum,  so  that  the  higher  APE,  the  higher  the  blue 
content  of  light.  A  spectroradiometer  is  required  for  obtaining 
this  input. 

The  model  takes  into  account  the  direct  normal  irradiance 
(DNI),  the  air  temperature  (Ta ir),  the  wind  speed  (WS)  and  the  APE. 
Coefficients  of  the  neural  network  are  obtained  from  outdoor 
monitorized  data. 


4.5.  Model  of  F.  Almonacid  et  al.  (2013) 

The  solar  research  group  of  Jaen  University  has  wide  experi¬ 
ence  on  the  application  of  artificial  neural  networks  in  the 
photovoltaic  field  [51-56],  Taking  into  account  this  experience 
and  from  a  similar  point  of  view  than  the  previously  commented 
model,  the  model  proposed  by  F.  Almonacid  et  al.  [57]  tries  to 
characterize  the  relation  between  atmospheric  parameters  and 
module  output  maximum  power  through  artificial  neural  net¬ 
works.  In  this  case,  a  feed-forward  neural  network  trained  with 
the  Levenberg-Marquardt  back-propagation  algorithm  was  used. 

The  model  takes  into  account  the  spectral  influences  through 
easily  measurable  parameters:  air  mass  (AM)  and  precipitable 
water  (PW).  This  way,  it  avoids  the  need  of  spectroradiometer. 
Inputs  of  the  model  are  the  direct  normal  irradiance  (DNI),  the  air 
temperature  (Tair),  the  wind  speed  (WS),  the  AM  and  the  PW. 
Coefficients  of  the  neural  network  are  obtained  from  outdoor 
monitorized  data. 


Table  1 

Comparison  of  characteristics  of  models  for  high  concentration  photovoltaic  cells:  required  inputs  and  main  advantages. 


Models: 

Equivalent  circuit 

Kinsey 

Dominguez 

Fernandez 

Sandia 

Required  inputs 

Concentration  ratio 

V 

V 

V 

V 

V 

Cell  temperature 

V 

V 

V 

V 

V 

Incident  spectrum 

Subcells  photocurrents 

Air  mass 

V 

V 

V 

V 

V 

Advantages 

Provides  the  whole  I-V  curve 

V 

V 

V 

Spectral  effects  quantified  with  high  accuracy 

V 

V 

V 

V 

Explicit  equations 

V 

V 

V 

V 

Algorithms  for  extracting  the  model  parameters  are  easy  to  implement 

V 

V 

V 

It  does  not  require  knowing  in  detail  the  incident  spectrum 

V 

V 

V 

It  does  not  require  indoor  characterization  of  the  solar  cell 

V 
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4.6.  Model  of  E.F.  Fernandez  et  al.  (2012) 

The  model  of  E.F.  Fernandez  et  al.  [35]  uses  air  mass  (AM)  in 
order  to  quantify  the  spectral  influences  on  HCPV  module  max¬ 
imum  power.  Therefore,  it  avoids  the  need  of  spectroradiometer  or 
isotype  cells.  The  aim  is  to  have  a  model  based  on  easily  measur¬ 
able  parameters.  The  main  equations  of  the  model  are  as  follows: 

Pmpp  =AiDNI  +A2Tair  for  AM<2  (17) 

PMpp  =  A}  DNI  +  A2Tair  T  A3AM  for  AM  >2  (18) 

These  equations  consider  that,  for  large  values  of  AM,  PMpp  is 
mainly  influenced  by  DNI,  Tair  and  AM,  while  for  small  values  of 
AM,  the  influence  of  AM  can  be  neglected  [58,59].  Coefficients  Alt 
A2  and  A3  are  obtained  by  regression  analysis  from  outdoor 
monitorized  data.  At  present,  the  model  has  been  only  evaluated 
at  the  south  of  Spain. 

4.7.  Model  of  Sandia  National  Laboratories  (2004) 

The  models  above  only  calculate  the  maximum  power  of  the 
module.  When  more  information  about  the  device  I-V  curve  is 
required,  an  alternative  is  to  use  the  model  of  Sandia  National 
Laboratories  [36,60-63],  which  calculates  five  points  of  the  I-V 
curve  as  shown  in  Fig.  4. 

Equations  of  the  model  are  formulated  as  a  function  of  the 
"effective  irradiance”,  which,  for  the  case  of  a  HCPV  module,  is 
defined  as  it  is  shown  in  Eq.  (12). 

The  model  equations  that  allow  calculating  the  maximum 


power  of  a  HCPV  module  are: 

/j(AM)  =  a0  +  a]  AM  +  a2AM2  +  a3AM3  +  a4AM4  (19) 

Bef  =  [DN1/1(AM)]/DNI*  (20) 

S  =  [mk(Tc  +  273.1 5)]/q  (21) 

Avmpp  =  Avmpp  0  T  mVnipp(  ^  ^ef )  (22) 

Wp  =  (C0Bef  +  CjB^f  )[/'mpp  +  aimpp(T'c-fc)]  (23) 


^  MPP  —  V  [pipp  +  C2NsSln(Bef )  +  C3Ns[<5ln(Bef)]2  + /*vmpp(  Tc_T*)  (24) 


PMPP  = /mppVmpp  (25) 

Eq.  (19)  approximates  the  spectral  correction  factor,  /,(AM), 
as  defined  in  the  standard  ASTM  E  973  [40],  by  a  fourth  order 
polynomial.  Eq.  (20)  calculates  the  effective  irradiance  (Ber).  that  is, 
the  irradiance  at  which  the  cells  actually  respond.  This  irradiance 
is  the  direct  normal  irradiance,  DNI,  corrected  with  the  spectral 
correction  factor,  /i(AM),  and  normalized  to  the  reference  irradi¬ 
ance,  DNI*.  Eq.  (21)  defines  the  S  parameter,  which  is  the  product 
of  the  effective  ideality  factor  of  the  multi-junction  cell  (m)  and 
the  thermal  voltage.  The  thermal  voltage  is  obtained  from  the 
Boltzmann  constant  (k),  the  electron  charge  (q)  and  the  cell 
temperature  (Tc).  Eq.  (22)  allows  the  determination  of  the  tem¬ 
perature  coefficient  /?vmPP-  This  coefficient  is  used  afterwards  for 
quantifying  the  effect  of  temperature  on  module  maximum  power 
point  voltage.  The  coefficient  is  expressed  as  a  linear  function  of 
the  effective  irradiance,  i.e.  it  is  allowed  to  vary  with  the  concen¬ 
tration  ratio.  Eqs.  (23)  and  (24)  calculate  the  module  maximum 
power  point  current  (/MPP)  and  voltage  (VMPP)  from  their  values  at 
reference  conditions  (/mpp*,  Vmpp*)-  The  reference  conditions  are 
defined  in  the  model  by  DNI*  Tc*  and  AM1.5.  aImpp  is  the 
temperature  coefficient  for  JMPP  and  Ns  is  the  number  of  cells  in 
series  for  the  module.  Finally,  Eq.  (25)  obtains  the  maximum 
power  of  the  HCPV  module. 

As  it  has  been  shown,  some  of  the  referenced  module  models 
use  the  mean  cell  temperature  of  the  module  (Tc)  as  input.  The 
direct  measurement  of  this  parameter  is  very  difficult,  due  to  the 
special  features  of  these  kinds  of  devices.  Because  of  this,  indirect 
methods  for  estimating  Tc  have  been  proposed  [42,64-69], 

Table  2  summarizes  the  main  characteristics  of  the  referenced 
models  for  HCPV  modules,  including  the  required  inputs  and 
advantages  of  each  one. 

Table  3  presents  a  quantitative  comparison  of  the  referenced 
models  for  HCPV  modules  characterization.  Four  models  have 
been  evaluated  at  the  Center  of  Advances  Studies  in  Energy  and 
Environment  of  Jaen  University  (ASTM  E  2527,  Almonacid, 
Fernandez  and  Sandia).  The  rest  of  the  models  have  been  eval¬ 
uated  by  their  respective  authors  and  results  can  be  found  in  the 
papers  that  describe  the  models.  With  respect  to  the  values  shown 
in  Table  3,  it  is  important  to  remark  that  a  direct  comparison  of  the 
provided  values  is  not  possible  because  the  different  authors  use 
different  measurement  and  evaluation  procedures.  Furthermore, 


Table  2 

Comparison  of  characteristics  of  models  for  high  concentration  photovoltaic  modules:  required  inputs  and  main  advantages. 


Models: 

ASTM  E  2527 

Peharz 

Steiner 

Rivera 

Almonacid 

Fernandez 

Sandia 

Required  inputs 

Direct  normal  irradiance 

V 

V 

V 

V 

V 

V 

V 

Ambient  temperature 

V 

V 

V 

V 

V 

Wind  speed 

V 

V 

V 

V 

Precipitable  water 

V 

V 

Cell  temperature 

V 

V 

Incident  spectrum 

V 

Subcells  photocurrents 

V 

V 

Average  photon  energy 

V 

Aerosol  optical  depth 

V 

Air  mass 

V 

V 

V 

Tracking  error 

V 

Advantages 

It  quantifies  spectral  effects 

V 

V 

V 

V 

V 

V 

Spectral  effects  quantified  with  high  accuracy 

V 

V 

V 

It  does  not  require  the  cell  temperature 

V 

V 

V 

V 

V 

It  does  not  require  spectroradiometer 

V 

V 

V 

V 

V 

It  does  not  require  indoor  characterization 

V 

V 

V 

V 

V 

It  provides  electrical  parameters  other  than  PMpp 

V 

V 

Easy  to  implement 

V 

V 

V 

V 

It  does  not  require  advanced  knowledge  on  artificial  neural  networks 

V 

V 

V 

V 

V 

Low  computational  effort 

V 

V 

V 

V 

V 

V 
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Table  3 

Errors  between  measured  and  calculated  maximum  power  of  the  models  for  HCPV 
modules  as  provided  by  different  authors. 


Model 

Parameter 

Value 

ASTM  E  2527 

RMSE 

4.6%* 

Peharz 

Absolute  RMSE 

1.3, 1.2, 1.6  and  0.6  W** 

Steiner 

NRMSE 

Between  2.5  and  3.2%*** 

Rivera 

MSE 

^***=f= 

Almonacid 

RMSE 

2.1%* 

Fernandez 

RMSE 

3.2%* 

Sandia 

RMSE 

3.4%* 

*  Mean  value  of  the  root  mean  square  error  for  two  HCPV  modules  measured 
during  2011  and  2012  at  the  Center  of  Advanced  Studies  in  Energy  and  Environ¬ 
ment  of  Jaen  University  (south  of  Spain). 

**  Values  of  the  absolute  root  mean  square  error  for  four  HCPV  modules 
measured  in  different  time  periods  from  2009  to  2012  at  the  Fraunhofer  1SE 
(Freiburg,  Germany)  [42],  The  nominal  power  of  the  modules  at  850  W/m2  of  direct 
normal  irradiance,  25  ”C  of  cell  temperature  under  the  incident  spectrum  AM  1.5d 
ASTM  G173-03  is  54.0,  50.1,  44.0  and  15.7  W  respectively. 

***  Range  of  the  obtained  normalized  root  mean  square  error  for  five  HCPV 
modules  measured  from  November  2011  to  March  2012  at  the  Fraunhofer  1SE 
(Freiburg,  Germany)  [44],  The  errors  are  normalized  to  the  nominal  power  of  each 
module. 

****  ya|ue  0f  the  mean  square  error  for  one  HCPV  module  measured  from 
November  2011  to  March  2012  at  Jaen  University  (south  of  Spain)  [48], 

the  reported  errors  are  expressed  by  different  parameters  in  each 
case.  However,  the  table  can  be  illustrative  of  the  order  of 
magnitude  of  the  errors  of  the  models. 

5.  Conclusions 

Nowadays,  we  can  find  in  the  scientific  literature  several  models 
for  the  electrical  characterization  of  HCPV  cells  and  modules.  These 
models  are  required  for  the  design,  monitoring  and  energy  prediction 
of  HCPV  systems.  A  review  of  the  existing  models  has  been  presented 
with  the  aim  of  helping  the  photovoltaic  professionals  and  research¬ 
ers  in  the  choice  of  the  most  suitable  model  for  each  specific 
application. 

Multijunction  solar  cell  models  calculate  the  relevant  electrical 
parameters  of  the  cell  as  a  function  of  irradiance,  spectrum  and  cell 
temperature.  In  order  to  evaluate  the  spectral  influences,  spectro- 
radiometer,  isotype  cells  or  a  solar  simulator  can  be  required. 
However,  we  can  find  other  simplified  models  that  require  less 
sophisticated  instruments  in  spite  of  a  little  loss  of  accuracy. 
Numerical  algorithms  are  often  needed  in  order  to  adjust  the  model 
parameters.  These  algorithms  are  not  always  easy  to  implement. 

HCPV  modules  models  calculate  the  relevant  electrical  para¬ 
meters  as  a  function  of  ambient  parameters  such  as  direct  normal 
irradiance,  air  temperature  or  wind  speed.  The  spectral  influences 
can  be  characterized  in  different  ways  (Z,  APE,  AM)  and  different 
equipments  are  required  to  obtain  these  spectral  parameters. 
Every  model  requires  outdoor  measurements  in  order  to  obtain 
the  model  coefficients  by  means  of  regression  analysis.  Some  of 
these  models  are  easy  to  implement  and  do  not  require  very 
sophisticated  instruments  to  be  applied  so  that  they  are  good 
candidates  to  be  used  in  the  industry,  in  particular  the  ASTM  E 
2527  [41],  the  model  developed  by  Fernandez  et  al.  [35]  and  the 
model  developed  by  Sandia  National  Laboratories  [36], 
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